Lecture 21 

Brahmagupta-Pell Equation 



Recall - For quadratic irrational x we defined 



xo = x= -^0 "j^ Co\d-BQ, d,Co,BoG 

ai = \xi\ 

B, + Vd 



Xi — 



Xi+l = 



Ci 
1 



Xi di 

Bi^i = aiCi — Bi 



Ci+i = 



Ci 



We showed that -Bj, Cj e Z, and that x has a purely periodic expansion if and 
only if a; > 1 and — 1 < x < 0. 

Corollary 72. Let dbe a positive integer, not a perfect square. Then the continued 
fraction of the number x = Vd + lVd\ is purely periodic. 

Proof. 

x = Vd+ [y/d] > 1 

X = Vd - [\/d\ satisfies - 1 <x <0 since [Vd\ < Vd < [Vd\ + 1 



Let's analyze thisx = Vd+ [Vd.\ a little more, x = ^^^+lV^^ and l\d- [Vd\'^, 
so we can take Cg = 1, and Bo = \^J. Want to see what happens for higher n - 
what Xn looks like. Let x = [ao, ai, . . . a^-i] be the continued fraction of x, r is 
chosen as smallest possible period. 

Claim: xq = x,xi,X2, ■ ■ ■ Xr-i are all distinct 

Proof. If Xq = Xi for some i < r — 1, then we'd have period i smaller than r ■ 

So Xn = Xo if and only if n is a multiple of r {xm = x„ if m = n mod r). We'll 
show that C„ = 1 if and only if n is a multiple of r, and C„ cannot be — 1. First, 
if n = fcr 



1 



Bkr + Vd 



c, 



kr 



= Xkr = Xn = Xo = Vd+ \ Vd\ 



Bkr — Ckr lVd\ — Vd{Ckr — 1) Only happens if Ckr = 1 (otherwise integer = 
irrational). Conversely, if C„ = 1 then x„ = Bn + Vd- 

We know Xn is also purely periodic [a„, a„+i, . . . a„+r-i]/ so 

Xn > 1 and — 1 < x^n < 

-1 < B„ - \/d < 
^ Bn <Vd< B„ + l 
^ B„ = [Vd\ 

which means that a;„ = Vd + lVd\ = xq, so that n is a miiltiple of r. 
Suppose C„ = — 1. Then a;„ = — B„ — Vd is purely periodic, so 

Xn> -Bn - Vd> I 

and 

-Kx:^ <0^ -1< -Bn + Vd<0 

which means that Bn > Vd and Bn < —Vd — 1 Vd < —Vd - 1, which is 
impossible. 

Note that ao = lx\ = [Vd + [v^JJ = lVd\ + lVd\ = 2lVd\. So continued 
fraction expansion of a; = Vd + lVd\ is 



[2[-\/cZJ , oi, . . . Ur-i] = [2[Vd\ , oi, . . . a^-i, 2[-\/dJ] 

Continued fraction expansion of Vd will look like that of Vd + lVd\ except with 
a different first digit [ [Vd,\ , ai , . . . a^— i , 2 [Vd\ ] . 

Note: We can run the {Bn,Cn) process for x ^ Vd ^ ^^,Co = 1, = 0, 
note that xi = -^rrj^ is the same for x = Vd and for x = Vd + [Vd\, so since 

Xn = '^"(j"^ is the same for these two x's as long as n > 1, and also because 
Xn = then Bn, Cn are the same forn > 1 whether we start with Vd or 

Vd + lVd\, so still true that Cn 7^ —1 and C„ = 1 if and only if n = kr. 



Theorem 73. IfdGN is not a perfect square, and {^} are the convergents to Vd, 

and C'n is the sequence of integers we defined for Xn (starting with xq = °~^/^ ), then 
pl-dql = i-ir+'Cn+,. 
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Proof. 



Xq 



dQn + Vd{Bn+iqn + g„_iC„_,_i) 



Xn+lPn + Pn- 



Xn+lQn + Qn-1 



Pn +Pn-1 



qn + qn-1 



\ C„_|_i J 
{Bn+lPn + Pn-lC„+l) + Vdpn 

(i?„+ig„ + g„_iC„+i) + Vdqn 

{Bn+lPn + Pn-lCn+l) + VdPn 



By comparing coefficients, we get that 

(Bn+iqn + qn-lCn+l)Pn = Pn 

{Bn+lPn + Pn-lCn+l)qn = dql 

Cn-1 {PnQn-1 " QnPn-l) = Pn - dq^ 
^ V ■' 

(_l)n-l 

pI - dql = {-ir+^Cn+i 



Corollary 74. Ifr is period of the continued fraction expansion of \fd, then p\^_x — 

Remark 2. If nr is even then we get a solution {pn,qn) of the P-B equation 
since pI^-i — dq^r-i = (—1)^^^ = 1, so we get infinitely many solutions since 
convergents are all distinct. 

Back to P-B equations - dy"^ = 1 with d e Z, want x,y e Z. li d < 0, then 
x'^ + \d\y'^ = 1, since x, y G Z, finite number of easily computed solutions. So, 
can assume d > 0. We showed last time that in fact, all solutions must come 
from continued fraction of \/d. 

More generally, (*) - dy"^ = N for N e Z. If {x, y) is a solution of (*), then 
so is (ia;, ±j/) for any choice of signs. Some trivial solutions for a; = or ?/ = 0, 
so look for nontrivial. Then we can assimie x,y > Q. These are called positive 
solutions. Also assume that {x, y) = 1. (If not, replace N with ^ i£ g = {x, y)). 
So only looking for positive, primitive (a;, y). 



Theorem 75. Let d € N, d ^ □, and let N e Z such that \N\ < Vd. Then 
any positive primitive solution {x, y) of x"^ — dy"^ = N has the property that | is a 

convergent to \/d. 

Proof. Suppose p is a positive real number such that is irrational and cr e M, 
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s, i e N such that -t^p = a and also that < cr < ,/p. 



Claim: 

1 

< 



2*2 



Proof of Claim. 

s s- ty/p 



t{s + t^p) 
_ - t^p 
~ t{s + t^) 
a 

Note that because — t^p = a > 0, s > t^, so s + > 2t^, so that 

□ 

Now, using the claim we see that f is a convergent to the continued fraction of 
y/p (by Problem 4 of PSet 9). 

U N > 0, just use a = N, p = d, (s, t) = {x, y) to show that | is a convergent 
to v^. If A'' < 0, rewrite x'^ - dy^ = N as j/2 - la;2 = - ^, then take a = - f . 
\N\ < \fd, so < (7 < ^ = and so | is a convergent to continued fraction 

Note that if the continued fraction of \fd = [ao, ai,... ], then continued fraction 
of ^ = [0, ao, ai, . . . ] means that convergents of ^ are just reciprocals of 

convergents of \/d. 

rn 1 1 1 gfc 

[0,ao,ai,---J = 1 = = — 

ao + ^ ^ Pk 

and so if ^ is a convergent to then ^ is a convergent to Vd ■ 

Theorem 76. Let d G N,d ^ H. All positive solutions to x^ — dy^ = ±lare of the 
form (x, y) = (p„, g„) a;/iere is convergent to \fd. Ifr is the period of the continued 
fraction of \fd, then 



4 



• Ifris even, — dy"^ = — 1 doesn't have any solutions, and all positive solutions 
ofx^ — dy"^ = 1 are given by x — pkr~i,y = qur-ifor = 1, 2, 3, 

• Ifris odd, then all positive solutions to x"^ — dy"^ = — 1 are given by taking x = 
Pkr-i,y = qkr-ifor k = 1,3,5, ... , and all positive solutions to — dy^ = 1 
are given by taking x = Pkr-iiV = Qkr-ifor k = 2,4,6, .. . 

Proof. If {x, y) is a positive solution to — dy'^ = ±1 then %cd{x, y) = 1 is forced. 
By theorem it must come from convergent to Vd, say . But we showed that 
- dq^^ = (— l)"+^C„+i. Also Cn+1 can't be —1, and can be 1 if and only if 
n + 1 is a multiple of r - ie., n = fcr — 1. So, Pfc^-i ~ '^llr-i — (— 1)*^*^ => if r even, 
can't be —1, and if r odd, can be ±1. ■ 

Remark 3. Suppose two positive solutions {xi,yi) and {x2, 2/2) are solutions of 
x'^ — dy"^ = 1, then xi < xi <s=^ y\ < yi. 

Proof, t/i < 2/2 => = 1 + dy\ < 1 + dy\ = X2 and xi, a;2 > so xi < X2. Same 
for other direction, which means that we can order the positive solutions ■ 

Theorem 77. If{xi , yi) is the least positive solution ofx^—dif = 1 where D ^ d G N, 
then all positive solutions are given by (a;„, t/„) where Xn + Vdyn = {xi + Vdyi)"'. 

Eg. For x"^ — = 1, (3, 2) is the smallest positive solution. Then (3 + 2^2)^ = 
17 + I2V2 (17, 12) is the next solution. 
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